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Abstract 

We present an updated discussion of K — > nil decays in a combined framework of chiral pertur- 
bation theory and Large- A^'c QCD, which assumes the dominance of a minimal narrow resonance 
structure in the invariant mass dependence of the II pair. The proposed picture reproduces very 
well, both the experimental K'^ vr^e^e" decay rate and the invariant e'^e^ mass spectrum. 
The predicted Br(_ft's —^ Tr'^e^e") is, within errors, consistent with the recently reported result 
from the NA48 collaboration. Predictions for the K -k h'^ fj,~ modes are also obtained. We 
find that the resulting interference between the direct and indirect CP-violation amplitudes in 
Kl — » 7r°e^e~ is constructive. 
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1 Introduction 



In the Standard Model, transitions like K ttI^I^, with I ~ e,fi, are governed by the interplay of 
weak non-leptonic and electromagnetic interactions. To lowest order in the electromagnetic coupling 
constant they are expected to proceed, dominantly, via one-photon exchange. This is certainly the 
case for the K"^ 7r='=Z+/" and Ks tt^H+I^ decays p. The transition 7r°7* 7r°/+Z~, via 

one virtual photon, is however forbidden by CP-invariance. It is then not obvious whether the physical 
decay Kl tt^^+Z" will still be dominated by the CP-suppressed 7*-virtual transition or whether a 
transition via two virtual photons, which is of higher order in the electromagnetic coupling but CP- 
allowed, may dominate The possibility of reaching branching ratios for the mode Kl 7r°e+e^ 
as small as 10""'^^ in the near future dedicated experiments of the NA48 collaboration at CERN, is a 
strong motivation for an update of the theoretical understanding of these modes. 

The CP-allowed transition K2 77^7*7* 7r"e+e~ has been extensively studied in the literature 
(see refs. |2| E] and references therein). We have nothing new to report on this mode. A recent 
estimate of a conservative upper bound for this transition gives a branching ratio ^ 

Br(/\L ^ 7r°e+e")|cpc < 3 x lO"^^ ^ (-^ ;L) 

There are two sources of CP-violation in the transition 7r°7* 7r°Z+Z~. The direct source 

is the one induced by the "electroweak penguin" -like diagrams which generate the effective local 
four-quark operators [H] 

Qii =4(si7^di) ^ ([^7^/^) and Q12 - 4 (5^7^^^) ^ ([^^7,^,^) (1.2) 

l—e.fL l—e.fL 

modulated by Wilson coefficients which have an imaginary part induced by the CP-violation phase 
of the flavour mixing matrix. The indirect source of CP-violation is the one induced by the K^- 
component of the Kl state which brings in the CP-violation parameter e. The problem in the 
indirect case is, therefore, reduced to the evaluation of the CP-conserving transition Tr'^e+e^. 
If the sizes of the two CP-violation sources are comparable, as phenomenological estimates seem to 
indicate |2E1IZ1E], the induced branching ratio becomes, of course, rather sensitive to the interference 
between the two direct and indirect amplitudes. Arguments in favor of a constructive interference have 
been recently suggested 

The analysis of K t:^* nl^l^ decays within the framework of chiral perturbation theory 
(xPT) was first made in refs. 0I2|- To lowest non trivial order in the chiral expansion, the corre- 
sponding decay amplitudes get contributions both from chiral one loop graphs, and from tree level 
contributions of local operators of 0{p'^). In fact, only two local operators of the 0{p'^) effective 
Lagrangian with AS* = 1 contribute to the amplitudes of these decays. With C^{x) the 3x3 flavour 
matrix current field 

C^{x) = -iF^U{x)^D^U{x) , (1.3) 

where U{x) is the matrix field which collects the Goldstone fields (tt's, K^s and rf), the relevant effective 
Lagrangian as written in ref. is 

Lfi-\x) - KdKTs gg {tr (A/:^/:^) - [wi iv{Q\C^C,) + W2 iv{QC^XC.)] F^^^ + h.c. . 

(1.4) 

Here is a covariant derivative which, in the presence of an external electromagnetic field source 
Afj, only, reduces to DfjU{x) = dfjJJ{x) — ieAfj^{x)[Q ,U {x)] ; F'^'^is the electromagnetic field strength 
tensor; Fq is the pion decay coupling constant (Fq ~ 87 MeV) in the chiral limit; Q the electric charge 
matrix; and A a short-hand notation for the SU{2>) Gell-Mann matrix (Ae — iXj)/2: 

/ 2/3 0\ /000\ 
0= -1/3 , A= . (1.5) 

\ -1/3 / \ 1 / 
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The overall constant gg is the dominant coupling of non-leptonic weak transitions with AS ~ 1 and 
AI = 1/2 to lowest order in the chiral expansion. The factorization of gg in the two couplings wi 
and W2 is, however, a convention. 

For the purposes of this paper, we shall rewrite the effective Lagrangian in Eq. 1)1. 4|1 in a more 
convenient way. Using the relations 



Q = Q-li 



QX = XQ = — -A and 
and inserting the current field decomposition 



Q = diag.(l,0,0), / = diag.(l,l,l) 



C^ix) = L^{x) ~ eF^A^{x)A{x) , 



where 



L^{:x) = -iF^U\x)d^U{x) 
in Eq. (|1.4(l . results in the Lagrangian 



and 



A(x) = f/t(x)[Q,C/(x)], 



(1.6) 

(1.7) 
(1.8) 



le 

3^ 



- wa) tr (XL^L,,) + 3 Wa tiiXL^QL^) | + h 



(1.9) 



The Qii and Q12 operators in Eq. (|1.2() arc proportional to the quark current density (slj^^cIl) and, 
therefore, their effective chiral realization can be directly obtained from the strong chiral Lagrangian 
[ [sLl^dL) ^ {^^)23 to 0{p)]. Using the equations of motion for the leptonic fields d^F^„ = el-f^J, 
and doing a partial integration in the action, it follows that the effect of the electrowcak penguin 
operators induces a contribution to the coupling constant w only; more precisely 



Wi - W2j 



47ra 



[Cn(M') + Ci2(M')] 



(1.10) 



where Cii(/i'^) and Ci2(/i^) are the Wilson coefficients of the Qn and Q12 operators. There is a 
resulting ^-scale dependence in the real part of the Wilson coefficient Cn + C12 due to an incomplete 
cancellation of the GIM~mechanism because, in the short-distance evaluation, the w-quark has not 
been integrated out. This /z-dependence should be canceled when doing the matching with the long- 
distance evaluation of the weak matrix elements of the other four-quark operators; in particular, with 
the contribution from the unfactorized pattern of the Q2 operator in the presence of electromagnetism. 
It is in principle possible, though not straightforward, to evaluate the w and W2 couplings within the 
framework of Large-iVc QCD, in much the same way as other low-energy constants have been recently 
determined (see e.g. ref. |H] and references therein). While awaiting the results of this program, we 
propose in this letter a more phenomenological approach. Here we shall discuss the determination of 
the couplings w and W2 using theoretical arguments inspired from Large- A^c considerations, combined 
with some of the experimental results which are already available at present. As we shall see, our 
conclusions have interesting implications for the CP-violating contribution to the Kl 7r"e+e^ 
mode. 



2 K ^ ttU Decays to 0(j9^) in the Chiral Expansion 

As discussed in ref. PP, at C'(p^) in the chiral expansion, besides the contributions from the Wi and W2 
terms in Eq. H1.4|l there also appears a tree level contribution to the 7r+e"'"e~ amplitude induced 

by the combination of the lowest 0{p'^) weak AS* = 1 Lagrangian (the first term in Eq. (|1.4)l ) with 
the Lg-coupling of the ©(p^) chiral Lagrangian which describes strong interactions in the presence of 
electromagnetism [Hj: 

= -ieL9F'^''(x)tr {Q D^U{x)D,U^x) + Q D^U{x)^ DM{x)} . (2.1) 
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In full generality, one can then predict the /v+ —>■ tt^I^I decay rates {I = e,^) as a function of the 
scale-invariant combination of coupling constants 

w+ = -3(471)2 _ ^ 3(^^ _ 4L9)] - - log — ^ , (2.2) 

where Wi, W2 and Lg are rcnormalizcd couplings at the scale v. The coupling constant Lg can be 
determined from the electromagnetic mean squared radius of the pion jEI]: Lg(Afp) = (6.9±0.7) x 10~^. 
The combination of constants W2 — 4L9 is in fact scale independent. To that order in the chiral 
expansion, the predicted decay rate r{K^ — > 7r+e+e~) as a function of w+ describes a parabola. The 
intersection of this parabola with the experimental decay rate obtained from the branching ratio 

Bt{K+ n+e+e-) = (2.88 ± 0.13) x 10"^ , (2.3) 

gives the two phenomcnological solutions (for a value of the overall constant gg = 3.3): 

w+ = 1.69 ±0.03 and vi^+ = -1.10 ± 0.03 . (2.4) 

Unfortunately, this twofold determination of the constant w+ does not help to predict the Ks — > 
7r°e+e~ decay rate. This is due to the fact that, to the same order in the chiral expansion, this 
transition amplitude brings in another scale-invariant combination of constants: 

1 1 Af 2 

w,, = --(47r)2[wi-W2]--log^. (2.5) 

The predicted decay rate r{Ks ^ 7r"e+e^) as a function of is also a parabola. From the recent 
result on this mode, reported by the NA48 collaboration at CERN jl2j : 

Br {Ks 7r°e+e-) = [5.8t2:3(stat.) ± 0.8(syst.)] x 10"^ , (2.6) 

one obtains the two solutions for Wg 

w, = 2.56t°:^° and w, = -1.90t°:^3 . (2.7) 

At the same 0{p'^) in the chiral expansion, the branching ratio for the ir^e^e'^ transition 

induced by CP-violation reads as follows 



Br (Kt 



(2.4 ±0.2) 



/ImAt 



CPV 

2 



(3.9 ±0.1) 



- Ws 



(3.1 ±0.2) 



ImA, fl 



10- 



X 10" 



(2.8) 



Here, the first term is the one induced by the direct source, the second one by the indirect source 
and the third one the interference term. With ^3] ImAt = (1.36 ± 0.12) x 10"'', the interference is 
constructive for the negative solution in Eq. (|2.7|l . 

The four solutions obtained in Eqs. (|2.4() and (|2.7(l . define four different straight lines in the plane 
of the coupling constants W2 — 4Lg and w (= wi — W2), as illustrated in Fig. 1 below. We next want 
to discuss which of these four solutions, if any, may be favored by theoretical arguments. 



3 Theoretical Considerations 

3.1 The Octet Dominance Hypothesis 

In ref. it was suggested that the couplings wi and W2 may satisfy the same symmetry properties 
as the chiral logarithms generated by the one loop calculation. This selects the octet channel in the 
transition amplitudes as the only possible channel and leads to the relation 

W2 = 4Lg , Octet Dominance Hypothesis (ODH) . (3-1) 
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Fig. 1 The four intersections in this figure define the possible values of the couplings which, at 0{p^) 
in the chiral expansion, are compatible with the experimental input of Eqs. \2.!^) and \2.f)]) . The 
couplings wi, W2, and Lg have been fixed at the v = A/p scale and correspond to the value gg ~ 3.3. 
The cross in this figure corresponds to the values in Eqs. J^H.2(]\) and \S.21]) discussed in the text. 

We now want to show how this hypothesis can in fact be justified within a simple dynamical 
framework of resonance dominance, rooted in Large~A^c QCD. For that, let us examine the field 
content of the Lagrangian in Eq. (|1.9() . For processes with at most one pion in the final state, it is 
sufficient to restrict A and to their minimum of one Goldstone field component: 



A = -i^— [$, Q] + . . . , and i,, = V2F^ 9^$ + • • • , 
■f"o 

with the result (using partial integration in the term proportional to ieg8W2) 



(3.2) 



/•AS=1 
•'-eff 



{x) 



^ KdK; §8 {^F^ tr (A5^$a^$) + 2F^A^ tr[A($Oa^a> - 9^ci>Q<f>)] 
iew2 a.F'^'^ tr[A($Qa^$ - a,,$Q$)] 

-h.c. (3.3) 



-ie^wF"" tr(A(9p$a^$ 



showing that the two-field content which in the term modulated by W2 couples to d^^F'^^ is exactly the 
same as the one which couples to the gauge field A^^ in the lowest 0{p'^) Lagrangian. As explained in 
ref. 0, the contribution to /\+ 7r+7(virtual) from this 0{p^) term, cancels with the one resulting 
from the combination of the first term in Eq. H3.3|) with the lowest order hadronic electromagnetic 
interaction, in the presence of mass terms for the Goldstone fields. This cancellation is expected 
because of the mismatch between the minimum number of powers of external momenta required by 
gauge invariance and the powers of momenta that the lowest order effective chiral Lagrangian can 
provide. As we shall next explain, it is the refiect of the dynamics of this cancellation which, to a first 
approximation, is also at the origin of the relation W2 = 4L9. 
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With two explicit Goldstone fields, the hadronic electromagnetic interaction in the presence of the 
term in Eq. (|2.1|1 reads as follows 



2L9 

-^0 



d.F^n tr(Q$a^$) 



(3.4) 



The net effect of the Lg-coupling is to provide the slope of an electromagnetic form factor to the 
charged Goldstone bosons. In momentum space this results in a change from the lowest order point 
like coupling to 

(3.5) 



^0 



In the minimal hadronic approximation (MHA) to Large- A'c QCD jl4j . the form factor in question is 
saturated by the lowest order pole i.e. the p(770) : 



Ml 



M2 



12 ' 



which implies Lg 



IL 

2M2 



(3.6) 



It is well known that this reproduces the observed slope rather well. 

By the same argument, the term proportional to W2 in Eq. H3.3|l provides the slope of the lowest 
order electroweak coupling of two Goldstone bosons: 



£ew(a;) = ~ie % 1/ ^K*, gg 2F^ ( - ) tr[A($Qa^<I> - d^^Qn 



(3.7) 



V2 V -0 

In momentum space this results in a change from the lowest order point like coupling to 

l^l-:r^Q^ (3.8) 



2Ff ■ 



Here, however, the underlying AS' = 1 form factor structure in the same MHA as applied to Lg, can 
have contributions both from the p and the i4'*(892) : 



1 



aMl 



ml. 



Af2 + Q2 



with a + /3 = 1 , 



(3.9) 



because at ^ the form factor is normalized to one by gauge invariance. This fixes the slope to 



W2 

2FI 



/3 



a 



Ml Ml, 



(3.10) 



If, furthermore, one assumes the chiral limit where Mp = Mk* , there follows then the ODH relation 
in Eq. H3.1|l : a result which, as can be seen in Fig. 1, favors the solution where both w+ and Wg are 
negative, and the interference term in Eq. 12. 8|) is then constructive. 



3.2 Beyond the C»(/) in xPT 



A rather detailed measurement of the e+e invariant mass spectrum in K 



decays was 



reported a few years ago in ref. |17j . The observed spectrum confirmed an earlier result |18| which 
had already claimed that a parameterization in terms of only w-|_ cannot accommodate both the rate 
and the spectrum of this decay mode. It is this observation which prompted the phenomcnological 
analyses reported in refs. [3|S]. Here, we want to show that it is possible to understand the observed 
spectrum within a simple MHA picture of Large- A^c QCD which goes beyond the 0{p'^) framework of 
xPT but, contrary to the proposals in refs. |3E1, it docs not enlarge the number of free parameters. 

We recall that, in full generality the AT^ 7r+e+e~ differential decay rate depends only on 
one form factor 4>{z): 



dV 
dz 



1244^)4^ l^'^'^-J 



1-4^ 
z 







(-4) 









(3.11) 
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where = z M'^ is the hivariant mass squared of the e^e pair, and 

r< M ^r* "^^^ rnt 



V2 



K 



K 



(3.12) 



The relation between 4>{z) and the form factor pfotted in Fig. 5 of ref. pTj . which we reproduce here 
in our Fig. 2 below for |/v(z)| , is 



\fv{z)\' = 



(3.13) 
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Fig. 2 Plot of the form factor \fv{z)\ defined by Eqs. UTTl]) and versus the invariant mass 

squared of the e~^e~ pair normalized to M^. The crosses are the experimental points of ref. \1 '1^ ; 
the dotted curve is the leading 0{p'^) prediction, using the positive solution for w+ in Eq. {2.4}) ; the 
continuous line is the fit to the improved form factor in Eq. I.V. below. 



The O(p^) form factor calculated in ref. ^ is 



with the chiral loop functions 



z) =\-w++(l)Kiz) + (l)^{z)\ , 



(3.14) 



^ / ^ 41 5 1 /4 ^ 

(bKiz) = --- + — + -[--l 

3 z 18 S \z 



arctan 




and (j)j,{z) = (j) [ z — ^ 
ml 



(3.15) 



The experimental form factor favors the positive solution in Eq. H2.4|l . but the predicted 0{p'^) form 
factor, the dotted curve in Fig. 2, lies well below the experimental points for z > 0.2. 

Following the ideas developed in the previous subsection, we propose a very simple generalization of 
the 0{p'^) form factor. We keep the lowest order chiral loop contribution as the leading manifestation 
of the Goldstone dynamics, but replace the local couplings W2 — 4Lg and w = wi — W2 in w_|- by the 
minimal resonance structure which can generate them in the z-channel. For W2 — 4L9 this amounts 
to the replacement: 
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2Fl / Ml Ml Ml, Ml 
W2-4L9^-f a-— ^bT-+/3^ ^ ' 



'^^(A/^-A/iz)(A4.-M|z) ' 



while for w it simply amounts to the modulating factor: 



M 



Af 2 - Af2.z 



(3.17) 



Notice that in the chiral limit where Mp = Mk*, F^^ Fq, and when z ^ 0, wc recover the usual 
0(p^) couplings with the ODH constraint W2 = 4L9. In our picture, the deviation from this constraint 
is due to explicit breaking, induced by the strange quark mass, and results in an effective 



2F^ I A/2 

w.-4L,.-^Ml-^). (3.18) 



More explicitly, the form factor we propose is 



_ Gs {{An) 



A/2 Ml - Ml, 



A/2 - Aiy (A/2 - A/^z) (Mf,, - A/|z) 

6 y M^ J 3 60 

where the first line incorporates the modifications in Eqs. H3.16|l and l|3.17|l . while the second line is 
the chiral loop contribution of ref. PP, renormalized a,t v = Mp, and where we have only retained 
the first two terms in the expansion of (j)K{z), while x(z) = (j)Tr{z) — (j)Tr{0)- With w and /3 left as 
free parameters, we make a least squared fit to the experimental points in Fig. 2. The result is the 
continuous curve shown in the same figure, which corresponds to a xLn = 13.0 for 18 degrees of 
freedom. The fitted values (using gg = 3.3 and = 92.4 MeV) are 

w = 0.045 ±0.003 and /3 = 2.8 ±0.1; (3.20) 

and therefore 

W2 -4Lg = -0.019 ± 0.003. (3.21) 

These are the values which correspond to the cross in Fig. 1 above. 

As a test we compute the K'^ — > Tr+e^e^ branching ratio, using the form factor in Eq. (|3.19l) with 
the fitted values for w and /?, with the result 

Bt{K+ -> TT+e+e") = (3.0 ± 1.1) x 10"^ , (3.22) 

in good agreement (as expected) with experiment result in Eq. 1)2.3(1 . 
The fitted value for w results in a negative value for w,, in Eq. 1)2.5(1 

Ws = -2.1±0.2, (3.23) 

which corresponds to the branching ratios 

Br {Ks 7r°e+e-) = (7.7 ± 1.0) x 10"^ , (3.24) 
Br {Ks ^ ^°e+e-) |>i65McV = (4.3 ± 0.6) x 10^^ . (3.25) 
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This is to be compared with the recent NA48 results in Eq. (|2.6|l and ^21 

Br {Ks -> TT°e+e-) |>i65McV = [3tJJ(stat.) ± 0.1(syst.)] x 10"^ . (3.26) 
The predicted branching ratios for the K tt IJ-^ fJ.^ modes are 
Br(A'+ ^ 7r+^+/^-) = (8.7±2.8) X 10"^ and Br (A's -> ttV+Z^") = (1-7 ± 0.2) x 10-^ (3.27) 
to be compared with 

Br (A'+ ^ n+^i+fi-) = (7.6 ± 2.1) x 10"** , ref. [TTI (3.28) 

Br {Ks ttV+A^") = [2.9ti:2(stat.) ± 0.2(syst.)] x 10"^ , ref. HH . (3.29) 

FinaUy, the resulting negative value for Wg in Eq. (|3.23|) . implies a constructive interference in 
Eq. H2.8|l with a predicted branching ratio 

Br {Kl 7r"e+e-) |cpv = (3.7 ± 0.4) x 10"" , (3.30) 

where we have used [El ImAt = (1.36 ±0.12) x lO"'^ and we have taken into account the effect of the 
modulating form factor in Eq. (j3.17() . 

4 Conclusions 

Earlier analyses oi K ^ tt e+e~ decays within the framework of xPT have been extended beyond the 
predictions of C'(p^), by replacing the local couplings which appear at that order by their underlying 
narrow resonance structure in the spirit of the MHA to Large- TVc QCD. The resulting modification of 
the C(p^) form factor is very simple and does not add new free parameters. It reproduces very well both 
the experimental decay rate and the invariant e~^e~ mass spectrum. The predicted IiT{Ks — > ir'^e'^e") 
and Br(Xs — > 7r°/^+/x~) arc, within errors, consistent with the recently reported result from the NA48 
collaboration. The predicted interference between the direct and indirect CP-violation amplitudes in 
Kl 7r°e^e~ is constructive, with an expected branching ratio (see Eq. (|3.3()|l ') within reach of a 
dedicated experiment. 
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